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Amorphous solids, as well as many other disordered materials, display remarkable universality 
in their low temperature acoustic properties. This universality is attributed to the attenuation of 
phonons by tunneling two- level systems (TLSs), facilitated by the interaction of the TLSs with 
the phonon field. Recently a Two-TLS model was introduced, where it was shown that inversion 
symmetric TLSs, i.e. TLSs having their two states related by local inversion symmetry, dictate 
phonon attenuation at low temperatures, and are responsible for the universal phenomena. The Two- 
TLS model suggests resolution to some long standing questions regarding the universal phenomena, 
such as the smallness and universality of the tunneling strength, and the energy scale of 3K below 
which universality is observed. A crucial part of the Two-TLS model lies in the magnitude and 
spatial dependence of effective TLS-TLS interactions involving inversion symmetric TLSs. Here we 
calculate numerically, using conjugate gradients method, the effective elastic TLS-TLS interactions 
for inversion symmetric and asymmetric TLSs, for pure and disordered solids, in two and three 
dimensions. We verify analytical estimates used in the Two-TLS model. We further show that 
the disorder induced power reduction of the spatial dependence of the interaction persists to short 
distances, not much larger than the interatomic distance. Our results provide strong support to 
the Two-TLS model, and contribute to the better understanding of the microscopic structure of 
amorphous solids and the low temperature properties of glasses. 



I. INTRODUCTION 

The existence of two-level tunneling defects as a 
generic property in amorphous systems was postulated 
four decades agop], ^ in an attempt to explain the re- 
markable universality in the low energy characteristics 
of amorphous solids as were found earlier by Zeller and 
Pohlg. In what is now known as the Standard Tunnel- 
ing Model (STM)[H, 0,01 the interaction of the tunneling 
two-level systems (TLSs) with the phonon field is given 
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solids in the universal regime, T < Tjj « 3K (notable 
exceptions are the phenomena of spectral diffusion [f| 
and the dipole gap arising at very low temperatures, 
T < 0.1K[7H9[). Reconciliation of the assumption of 
non-interacting TLSs with the presence of TLS-phonon 
interaction is explained by the low density of TLSs, and 
consequent smallness of the typical interaction. 

Recently a new model was suggested, which is based 
on the different properties of the interaction of inver- 
sion symmetric and inversion asymmetric TLSs with the 
phonon field [13]. A TLS in which the two states are 
related to each other by inversion with respect to a mid- 
point, e.g. 180° flip of CN impurity in a KBr lattice, 
does not interact with the phonon field, but only with 
the second derivative of the phonon amplitude @ 



Here S? = ±1/2 denote the two states of the TLS at 
site j, u a P = du a /dr^ is the phonon field, where u is 
the phonon amplitude, a, f3 = x, y, z, and 7 is the TLS- 
phonon interaction parameter. The interaction of the 
phonon field with the tunneling amplitude of the TLS is 
small, and therefore neglected jl|, Q. 

The TLS-phonon interaction as given in Eq.(pQ) gives 
rise to an effective TLS-TLS interaction mediated by the 
phonon field [5|, [6j, which takes the form 
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where Uij is effectively random interaction with 1/r 3 de- 
pendence for distance r ^> clq. Here ao denotes the typ- 
ical interatomic distance. The STM, however, assumes 
that the TLSs are non-interacting, an assumption that is 
sufficient to explain many of the properties of amorphous 
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where we denote by r the symmetric TLSs, to differ from 
the asymmetric TLSs denoted by S. This form of the 
interaction leads, in three dimensions, to a long distance 
dependence of 1/r 4 and 1/r 5 for the acoustic mediated 
S TLS-r TLS (S-r) effective interactions and r TLS -r 
TLS (t-t) effective interactions respectively [6]. In two 
dimensions similar analysis to the one carried in Ref. [6[ 
leads to 1/r 2 , 1/r 3 , 1/r 4 spatial dependence of the S TLS 
-S TLS (S-S), S-r, and t-t interactions respectively. 

Where the above is correct for 2 impurities in an oth- 
erwise pure lattice, the presence of disorder introduces 
deviations from inversion symmetry, and a finite, albeit 
small interaction between a r-TLS and the phonon field. 
In Ref.[10| it was argued that for systems with strong 
disorder the TLS-phonon interaction can be written as 
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The small dimensionless parameter of the theory [10] is 
defined by g = 7 w /7s- This parameter is proportional 
to the deviations from inversion symmetry, and therefore 
to the ratio between the strain and interatomic lattice 
spacing, i.e. in strongly disordered systems g w 0.01 — 
0.03. The resulting effective elastic TLS-TLS interaction 
takes the formfioj 



the notion that each CN impurity constitutes a two-level 
system. Thus, the TLS-DOS in disordered lattices is not 
small, the relevant TLSs are r-TLSs, between which the 
interaction is small and can therefore be neglected, and 
the smallness of the DOS of detected TLSs is related to 
the fact that only TLSs with appreciable tunneling am- 
plitude can be detected. Note that in KBr:CN the TLS is 
formed by each CN low energy state and the state related 
to it by a CN flippy 

and is thus a r-TLS, where 
all other single particle states of a given CN impurity are 
much higher in energy, as they constitute an S excitation 
with respect to the CN low energy state [Io|. 
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where in three dimensions (two dimensions) all interac- 
tions decay as 1/r 3 (1/r 2 ) at distances r ^> ao, and their 
typical values at near neighbor distance are related by 

ur^gut^g 2 ul s . (6) 

It is further argued [lOj that the above long distance 
power law dependence of the interaction persists down to 
a short distance cutoff, which is not much larger than the 
interatomic spacing, for lack of another relevant length 
scale in the system. The form of the Hamiltonian (|5]), 
(|6]), i.e. the strength and spatial dependence of the three 
TLS-TLS interactions, and its subsequent analysis have 
led to a significant advance in the understanding of the 
low temperature universality in glasses [Toj. The univer- 
sality and smallness of phonon attenuation [ill [l2| were 
shown to be a consequence of the generic characteristics 
of the r TLSs, and the appearance of universality below 
Tjj w 3K was attributed to the gapping of the S-TLSs 
below the energy scale of Uq t « gU$ s w gTc ~ 3K, 
where Tq is the glass temperature. The enhancement of 
the single particle S-TLS DOS at higher temperatures 
results in their domination of acoustic attenuation, and 
the limiting of the universal characteristics of amorphous 
and disordered solids to T < Tjj. 

The Two-TLS model further derives, and quantifies, 
some of the central assumptions of the STM. The DOS 
of the relevant TLSs at low temperatures, the r-TLSs, 
is found to be finite at very low energies and rather ho- 
mogeneous for energies smaller than 3K. The assumption 
of noninteracting TLSs is supported by the fact that at 
low temperatures the S-TLSs are frozen, and the typical 
r — r interaction Uq t « 0.1K is smaller by a factor of 
g in comparison to the typical r-TLS energy, which dic- 
tates the typical energy scale of the universal phenomena. 
Thus, within the Two-TLS model TLS-TLS interactions 
can be neglected except at very low temperatures. How- 
ever, this possibility to neglect the interactions does not 
depend on the small concentration of TLSs, but on the 
form of the Hamiltonian ([5]), (|6]). This is in line with 
the existence of universal phenomena in disordered lat- 
tices such as KBr:CN, where there are strong evidences to 



The Two-TLS model, and specifically the Hamiltonian 
(0, was derived microscopically for strongly disordered 
lattices. However, the strong evidence for the equiva- 
lence of the phenomenon of the low temperature univer- 
sality in disordered lattices and amorphous solids (as well 
as other disordered systems, e.g. disordered polymers, 
metallic glasses, quasicrvstals) [l2i [liSl [l7|, strongly sug- 
gests that the mechanism leading to universality in dis- 
ordered lattices and amorphous solids is the same. Thus, 
validation of the Two-TLS model is significant both for 
the resolution of the long standing problem of univer- 
sality in disordered and amorphous solids, and for the 
advance of our understanding of the microscopic struc- 
ture of amorphous solids. The two crucial points of the 
Two-TLS model are the form of the Hamiltonian (|5]), 
(j6j), and the resulting structure of the DOS of the r 
and S TLSs. With regard to the latter, the results in 
Ref. [Toj| are supported by a comprehensive numerical cal- 
culation given in Ref. [18]. With regard to the former, a 
significant step was made in Ref. jl5j , where 7 w ,7s were 
explicitly calculated for KBr:CN. The obtained ratio of 
7w/7s ~ 0.02 and value 7 W w O.leV are both in agree- 
ment with theory [1Q] and with the experimental value 
for the coupling constant for the relevant TLSs at low 
temperatures. [l9[ l20j . These results support both the 
Two-TLS model and the notion that it is indeed CN 
flips that constitute the TLSs for KBr:CN. However, the 
ab-initio and DFT calculations [15| required the use of 
very small samples, not allowing the study of the effective 
TLS-TLS interactions and their distance dependence. 



In this paper we use the conjugate gradients method to 
study the magnitude and spatial dependence of the elas- 
tic impurity-impurity interactions in KBr:CN. Since CN 
impurities possess also charge, the elastic energy is de- 
rived by calculating first the total interaction energy and 
subtracting from it the direct electric TLS-TLS interac- 
tion. Calculations are performed for S-S, S-r, and r-r 
interactions in two and three dimensions, in both pure 
and disordered lattices, as function of distance between 
the TLSs. In Sec. [II] we present the calculation setup. 
Our results are presented in Sec. IIII1 and conclusions in 
Sec. El 
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A aa (kJ/mol) 


a aa (1/A) 


B aa (A 6 kJ/mol) 


KK 


158100 


2.985 


-1464 


CC 


259000 


3.600 


-2110 


NN 


205020 


3.600 


-1803 


BrBr 


429600 


2.985 


-12410 



TABLE I: Interatomic potential parameters (taken from Ref. 
[211] V Cross-interaction parameters were calculated from the 



combining rules A a p = (A aa Ap 

B a — (BaaBpp) 1 / 2 . 
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II. THE MODEL 

All our calculations are carried for CN impurities in 
KBr lattice. For some even size N, we arrange initially 
NxN K + and Br' ions in the two dimensional grid (or 
NxNxN ions in 3D grid) having the distance r = 3. 2 735 A 
between K + and Br~ ions [2l| . Then we replace some of 
the Br~ ions with CN" ions. The if + and Br ions are 
assumed to carry +1 and -1 charges respectively, while 
the charge of the CN' ion is represented by fractional 
charges qc = —1.28 and qAr = —1.37 placed on the 
carbon and nitrogen atoms, and the additional charge 
Reenter = +1.65 is placed at the center of mass [22]. The 
C-N bond length is fixed at 1.1 7 A, while the distances of 
the carbon and nitrogen atoms from the center of mass 
are fixed at 0.63A and 0.54A respectively [22|. Inter- 
atomic potential is calculated by the formula: 
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The interatomic potential parameters A a p, a a p and 
B a p are taken from Ref. [2l| and shown in Table HI and 
K = 1389.35 AkJ/mol. 

Non-linear conjugate gradients method (Fletcher- 
Reeves) is used to find the closest local minimum of the 
system potential energy U. Periodic boundary conditions 
are used to simulate the infinite crystal. 

In order to study the spatial dependence of the CN-CN 
interactions in the KBr:CN crystal we calculate the po- 
tential energy U r (as described above) for KBr:CN grids, 
containing only two CN ions. The CNs are placed at some 
distance R from each other in some direction [x, y, z\. For 
example, if the first CN ion is placed at coordinate [0,0,0] 
in the grid (3D array), R = 3 and the direction is [2,2,2], 
then the second CN is placed at [R • x, R • y, R • z] = 
[6, 6, 6]. For a given distance R and direction [x, y, z] we 
calculate the energies for several different orientations of 
CN ions in the grid and we use these energies to eval- 
uate the S-S, S-r and r-r interactions, as explained in 
detail below. Because the non-linear conjugate gradients 
method is very time consuming for large samples, we use 
both two dimensional and three dimensional grids. The 
three dimensional grids relate to the relevant experimen- 
tal systems, whereas the two dimensional grids allow us 



to study spatial dependence of CN-CN interactions at 
longer distances. 
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FIG. 1: S-S, S-r and r-r interactions calculation in 2D grid 
with two CN ions placed at some distance r from each other. 
Each pair of CN ions in the figure represents the energy Ur of 
the low energy state, where the ions are placed in the direction 
shown in the figure. The interaction energy U for each one 
of the interactions is calculated by U=(U r (l)-U r (2))-(U r (3)- 
U r (4)). 

We now describe the algorithm used to calculate the 
S-S, S-r and r-r interactions. We start by finding the 
low energy states for a single CN impurity in an other- 
wise pure KBr:CN lattice. In two dimensions these low 
energy states are four- fold degenerate, along the square 
diagonals. In three dimensions these states are eight-fold 
degenerate along the in-space diagonals. We emphasize 
that although for a single CN impurity the ground state 
degeneracy is lifted by tunneling, for strong disordered 
systems (CN concentration 0.2 < x < 0.7) the CN-CN in- 
teractions are much larger than the tunneling amplitude, 
and therefore it is the interaction with neighboring CN 
impurities that breaks the single impurity ground state 
degeneracy. The treatment of tunneling can then be de- 
ferred until the bias energies are determined [To|- Since 
in this paper we are concerned only with the energy bias 
of the different orientational states resulting from inter- 
actions with other CN impurities, tunneling is neglected 
altogether. This is also in line with our chosen calcula- 
tion method, which does not allow tunneling over large 
barriers. 

Let us now describe the calculation of the interaction, 
starting with two dimensions, see also Fig. [TJ For the 
S-S interaction at some distance r we calculate first the 
elastic energy U® las (l). Both CN ions are placed initially 
at an angle ^1,2 = tt/4 (angles are calculated between the 
direction of the N-C vector and the x-axis). The whole 
system is then relaxed. The CN-CN interaction results 
in a small shift in the CN positions ( < 00) and an- 
gles ( <C 7r). The resulting total energy of the system, 
U* ot (l) and positions of the CN molecules are recorded. 
The direct electric energy between the two CN molecules 
U^ lec (l) is calculated from the relaxed positions of the 
CN molecules, and U* las (l) = U* ot (l)- U^ lec (l) is de- 
duced. We then follow the same procedure to calculate 
U® las (2) with first CN ion having 0i « 7r/4 and the sec- 
ond ion is rotated with 02 ~ 37r/4, U® las (3) with first 



CN ion rotated having 0i « 37r/4 and the second ion 
with 2 ~ tt/4 and U® las (4) having both CN rotated with 
0i,2 ~ 37r/4. Finally, the S-S interaction energy is calcu- 
lated by: U elas (S-S) = (U* las (l) - U* las (2)) - (U* las (3) - 
U® las (4)). S-r and r-r interactions energies are calcu- 
lated in a similar way. The difference in S-r and r-r 
interactions calculation is that U® las (3) and U® las (4) are 
calculated while the first CN ion instead of being rotated 
is flipped, i.e. 0i ~ 57r/4. Additionally, for the r-r calcu- 
lation the second CN in U^ las (2) and U^ las (4) is flipped: 
02 ~ 57r/4, see Fig. [U 



Direction 


1st CN 
position 


2nd CN 
position 


3rd CN 
position 


of 
3rd CN 


of 
3rd CN 


[1,1,0] 
[2,0,0] 


[0,0,0] 
[0,0,0] 


[r,r, 0] 
[2 • r, 0, 0] 


[r, r, 2] 
[2-r + l,3,0] 


tt/4 
3tt/4 


0.955 
0.955 



TABLE IV: Positions of three CN ions and the angles of third 
CN ion in the three dimensional grids, for directions [1, 1,0] 
and [2,0,0]. 



III. RESULTS 



Interaction 


Angle 


U r l 


U r 2 


U r 3 


U r 4 




0(CAT1) 


tt/4 


3tt/4 


tt/4 


3tt/4 


S-S 


0(CN1) 


0.955 


0.955 


0.955 


0.955 




0(CAT2) 


tt/4 


tt/4 


3tt/4 


3tt/4 




6(CN2) 


0.955 


0.955 


0.955 


0.955 




0(CiVl) 


tt/4 


3tt/4 


tt/4 


3tt/4 




0(CN1) 


0.955 


0.955 


0.955 


0.955 


S-r 












0(CAT2) 


tt/4 


tt/4 


5tt/4 


5tt/4 




6(CN2) 


0.955 


0.955 


2.186 


2.186 



TABLE II: Approximate values of and angles of two CN 
ions in the three dimensional grids, used in the calculation of 
U r (i) for S-S and S-r interactions. 



In three dimensions the calculations are carried in a 
similar way. Since the low energy states of a single CN 
impurity are eight-fold degenerate, each CN impurity, 
whereas it can assume only a single flip, can rotate in 
six different directions. We thus specify for each calcu- 
lation which of the allowed rotations is performed, see 
Table HI 

As mentioned above, the spatial dependence of the S-r 
and r-r elastic interactions is modified in the presence of 
disorder. We therefore repeat the above calculations for 
CN-CN interactions in the presence of a third and fourth 
impurities. As the purpose of these added impurities is 
to introduce strain disorder, they are placed in a given 
position and orientation, and are then relaxed with the 
whole system for each configuration of the two "original" 
CNs between which the interaction is calculated. 



No 


1st CN 


2nd CN 


3rd CN 


4th CN 


of 


of 


position 


position 


position 


position 


3rd CN 


4th CN 


1 


[0,2] 


[2-r,2] 


[2-r + l,3] 




tt/4 




2 


[0,2] 


[2t,2] 


[2t + 1,5] 


[2-r,0] 


3tt/4 


tt/4 



TABLE III: Positions of four CN ions and the angles of third 
and fourth CN ions in the two dimensional grids, for direction 
[2, 0]. First row corresponds to the data denote "3 CNs 1" in 
Figs. [4] [5] Second row with the first three and with all four 
impurities correspond to the data denoted "3 CNs 2" and "4 
CNs 2" in the above figures respectively. 



The S-S, S-r, and r-r elastic interaction energies 
[U elas (S-S), U elas (S-r) and U elas (r-r)] for two CN im- 
purities in an otherwise pure lattice calculated for two 
dimensional grids of size 50X50 are plotted in Fig. [2] 
and Fig. [3] for directions [2,0] and [3,1], respectively. 
The 1 unit of distance in the plot is chosen to be the 
distance between two closest ions in the grid, which is 
about 3.2735A. Calculations are performed up to maxi- 
mal distance of 20 ions in diagonal, as at larger distances 
boundary condition effects become significant. Data at 
the shortest distances (<~ ao) is disregarded, as it shows 
irregularity dictated by the properties of the lattice at 
near neighbor distance. 

Since our available distance range is limited, lattice 
discreteness is expected to play a role. Indeed, the ex- 
act values of the magnitude of the elastic interactions at 
short distance and their distance dependence depends on 
details such as the direction between the two CN impu- 
rities, and the detailed position and orientation of the 
additional impurities, see below. Still, our numerical re- 
sults strongly support both the typical magnitude of the 
various elastic interactions (S-S, S-r, r-r) at short dis- 
tance, and their spatial dependence, as discussed above. 

The magnitude of the elastic interactions we find nu- 
merically is in agreement with experiments, and with 
their interpretation by the Two-TLS model. The elas- 
tic S-S interaction at short distance is of the order of 
100K, which is the same order of the glass temperature 
in KBr:CN. The elastic S-r interaction is of the order 
of 3K, which is g times smaller, and is characteristic of 
the temperature below which universality is observed flQ| . 
The elastic r-r interaction is of the order of 0.1K (another 
factor of g smaller), in agreement with the theory of the 
Two-TLS model, and with low energy experiments re- 
porting interaction phenomena at this energy scale [H, Q. 
The spatial dependence is also in good agreement with 
theory @, as the S-S, S-r, and r-r elastic interactions be- 
have as l/r a , with a « 2, 3, 4 respectively. We note that 
because the magnitude of the elastic r-r interactions at 
distances larger than a few times ao is not much larger 
than our numerical error, reliable results for the r-r in- 
teraction are available only for two dimensional samples 
in the (2,0) direction. 

The power law reduction of the r dependence of the 
elastic S-r interaction upon the introduction of disorder 
(from 4 to 3 in three dimensions, from 3 to 2 in two di- 
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FIG. 2: Spatial dependence of the elastic S-S, S-r and r-r 
interactions of 2 CNs placed in the direction [2,0] in 2D grid of 
size 50x50. Here and in subsequent figures energies are given 
in Kelvin units, l(KJ/mol)//c B = 120.274if. Solid line are 
best linear fits, where slopes denote the power law distance 
dependence of the interaction. 
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FIG. 3: Spatial dependence of the elastic S-S and S-r inter- 
actions of 2 CNs placed in the direction [3,1] in 2D grid of 
size 50x50. 



mensions) is a crucial aspect of the Two-TLS modelfioj. 
In order to test this result numerically, we introduce ad- 
ditional, third and fourth, impurities to the sample, see 
Table IIIII We then repeat the calculations above for the 
original two CN impurities. The additional impurities 
introduce strain disorder at the location of the measured 
TLSs. The results for the spatial dependence of the elas- 
tic S-r interaction with two, three and four CN impurities 
in 2D grid in direction [2,0] are shown in Fig. HI Despite 
the large variance in the calculated power laws for the 
different disorder configurations, a clear reduction of the 
power law is obtained, consistent with the theoretically 
predicted reduction from 3 to 2 [10]. Our results here also 
support the theoretical prediction that disorder does not 
affect the typical magnitude of the elastic S-r interaction 
at short distance. In Fig. [5] the overall results of S-S, S-r 
and r-r elastic interactions in the 2D grids are plotted. 



The factor of g ~ 0.03 separating the magnitudes of the 
S-S, S-r, and r-r interaction, the insignificance of disor- 
der for the typical magnitude of these interactions at the 
shortest distances, and at the same time the effect of dis- 
order on the spatial dependence of the S-r interactions, 
are all clearly demonstrated. A central result of this pa- 
per is the finding that the power law dependence of the 
various interactions, and specifically the disordered in- 
duced spatial dependence of the S-r interactions (~ 1/r 2 
in two dimensions, ~ 1/r 3 in three dimensions, see be- 
low), persists to very short distances, not much larger 
than the interatomic spacing. 
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FIG. 4: Spatial dependence of the elastic S-r interaction 
with two, three and four CNs placed in the direction [2,0] 
in 2D grid of size 50x50. Data for 2 CNs corresponds to 
S-r interaction in a pure KBr lattice. Data for 3 and 4 CNs 
corresponds to S-r interaction (calculated for impurities 1 and 
2) in a disordered lattice (introduced by impurities 3 and 4). 
Positions and orientations of third and forth impurities are 
given in Table Hill In Table Hill row 1 corresponds to the data 
denoted here "3 CNs 1" and row 2 to the data in "3 CNs 2" 
and "4 CNs 2". 
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FIG. 5: Spatial dependence of the elastic S-S, S-r and r-r 
interactions in 2D grids of size 50x50. 

Similar calculation, for the elastic S-S interaction in 
a pure lattice, and for the elastic S-r interactions in a 
pure lattice and with disorder strain induced by a third 
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FIG. 6: Spatial dependence of the elastic S-S and S-r interac- 
tions in 3D grids of size 20x20x20. Positions and orientations 
of impurities are given in Table IIVI 

impurity, were performed for three dimensional samples 
of size 20x20x20. Results for directions [1,1,0] and [2,0,0] 
with impurities in the positions specified in Table IIVI are 
shown in Fig. [6j For the elastic S-r interaction in an 
otherwise pure lattice we find for the spatial dependence 
a good agreement with the 1/r 4 behavior predicted by 
theory. For the power law of the spatial dependence of 
both the elastic S-S interaction and the elastic S-r inter- 
action with strain disorder we obtain values of 2.1 — 3.01. 
We attribute this large variance in exponents to the sig- 
nificance of lattice discreteness and boundary conditions, 
consequences of the small lattice size. Still, the reduction 
of the exponent of the spatial dependence of the elas- 
tic S-r interaction upon the introduction of disorder is 
clearly shown. With regard to the typical magnitude of 
the interactions, our results clearly reproduce the correct 
magnitude of order 100K for the elastic S-S interaction 
at short distances, and the much smaller, by factor of 



g w 0.03, magnitude of the elastic S-r interaction, with 
and without disorder. 



IV. CONCLUSIONS 

Using conjugate gradients method we have calculated 
the elastic interactions between rotations and flips of CN 
impurities in a KBr lattice. CN flips constitute two states 
related to each other by inversion symmetry, defined as 
r TLSs. CN rotations constitute asymmetric (S) TLSs. 
For two CN impurities in an otherwise pure lattice we 
find that the elastic S-S, S-r, and r-r interactions dif- 
fer both in their magnitude at short distances, and in 
their spatial power law dependence. Introducing disor- 
der does not change the relative magnitude of the in- 
teractions at short distances, but affects the spatial de- 
pendence of interactions involving a r-TLS, as inversion 
symmetry is broken. Our results are in agreement with 
theoretical predictions both for the magnitude and the 
power law dependence of the interactions. For the lat- 
ter, we find that the power law dependence predicted 
for large distances persists down to very short distances, 
not much larger than the interatomic distance. Since the 
magnitude and spatial dependence of the elastic S-r in- 
teractions is a crucial element in the Two-TLS model for 
the low temperature universality in disordered and amor- 
phous solids, our results contribute to the validation of 
the above model, and to the advance of our understand- 
ing of the microscopic structure of amorphous solids and 
of the low temperature properties of glasses. 
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